The number of real roots 
of a bivariate polynomial on a line 

Martin Avendano * 
February 2, 2008 

Abstract 

We prove that a polynomial / G R[a;, y] with t non-zero terms, restricted to a real line 
y = ax + b, either has at most 6t — 4 zeroes or vanishes over the whole line. As a consequence, 
we derive an alternative algorithm to decide whether a linear polynomial divides a sparse poly- 
nomial / G K[x,y] with t terms in [log{H{f)H{a)H{b))[K : Q] log(deg(/))t]°'^' bit operations, 
where if is a real number field. 



1 Introduction 

The famous Descartes' Rule of Signs, 1641, establishes that the number of positive real roots of a 
polynomial / G M[a;], counted with multiplicities, is bounded by the number of changes of signs in 
its ordered vector of coefficients, disregarding the zeros. As a direct consequence, the number of 
different real roots of / is bounded by 2t — 1, where t is its number of non-zero terms (here all roots 
are counted with multiplicities, except which is counted at most once). 

There are not yet natural generalizations of Descartes' Rule of Signs for the multivariate setting, 
but a lot of work has been and is being done for estimating the number of real isolated or non- 
degenerate roots (that is where the Jacobian does not vanish, condition that implies that the root 
is isolated) of multivariate square systems of real polynomials in the positive orthant, in terms of 
the number of variables and the number of non-zero terms that the system involves. 

The main result in that direction is due to A. Khovanskii 4 . A simple version of it implies that a 
square system of n real polynomial equations in n indeterminates, which involves in total t non-zero 
terms has at most {n + 1)*2**^*~^^/^ non-degenerate real roots in the positive orthant. Improvements 
of Khovanskii's bound have afterwards been obtained by D. Perrucci [TT] and T.Y. Li, J.M. Rojas 
and X. Wang [3], but for general systems the exponential dependence on the number of non-zero 
terms t can not be avoided yet. 



* Departamento de Matematica, Facultad de Cicncias Exactas y Naturales, Universidad de Buenos Aires, 
1428 Buenos Aires, Argentina. Research supported by grants CONICET PIP 2461/01 and UBACYT X-112. 
mavendarOdm . uba . ar 



1 



In [S], T.Y. Li, M. Rojas and X. Wang studied particular cases of bivariate square systems and 
showed that the number of common isolated or non-degenerate roots of a trinomial and a polynomial 
with at most t non-zero terms, t ^ 3, is bounded by 2* — 2. 

Furthermore, Kushnirenko's Conjecture, formulated in the mid-1970' (which says that a square 
system of n real polynomial equations in n indeterminates such that the A:-th polynomial has tk 
non-zero terms should have at most (ii — 1) • • ■ (t„ — 1) non-degenerate roots in the positive orthant) 
turned out to be false, by the counter-example provided by B. Haas in 2002 for a system of two 
trinomials in two variables [8J. 

The main result of this article is a refinement of the previous result for the particular case when the 
trinomial is a linear polynomial. Without loss of generality we can assume the linear polynomial is 
of the form y — ax — h and wc thus study the possible number of real roots of a bivariate polynomial 
on a line y = ax -\-b: 

Theorem 1.1 Let f = X]'=i o,ix"'y^^ G R[x, y] be a polynomial with at most t non-zero terms, and 
let a, 6 G M. Set g(x) — f{x, ax + b). Then either g = or g has at most 6t — 4 real roots, counted 
with multiplicities except for the possible roots and —b/a that are counted at most once. 

To our knowledge, this is the first time a non-exponential bound is achieved, even for systems of 
very particular shape like this one. The tools we use are completely elementary, and we are now 
studying the possibility of extending the results for more general systems. 

As a consequence of our result we derive an alternative algorithm for checking if a given linear form 
y — ax — b divides a polynomial / in in K[x, y], where K is a real number field. The number of bit 
operations performed by the algorithm is polynomial in the degree [K : Q] of the field extension, 
in the number t of non-zero terms of /, in the logarithm of the degree of / and in the logarithmic 
height of a, b and /. 

The first algorithm for this purpose can be deduced from a more general result by E. Kaltofen 
and P. Koiran |10| . They showed a polynomial-time algorithm for computing all linear factors of a 
sparse bivariate polynomial. This result has been further generalized in [12j and [13] to an algorithm 
that computes all the small degree factors of bi- and multi-variate sparse polynomials. All these 
algorithms use a version of the "gap theorem" introduced by F. Cucker, P. Koiran and S. Smale [SI. 
Instead of it, we reduce the problem to the univariate case by considering specializations f{x,x") 
for small values of n. 

2 Proof of Theorem 11.11 

Definition 2.1 Let f ~ X]f=o ^ R[x] be a non-zero polynomial. We note by V{f ) the number 
of changes of signs in the ordered vector (od, . . . , ag) of the coefficients of f , disregarding the zeroes. 
We also set ¥(0) = -2. 

Next result is our crucial ingredient in the proof of Theorem 11.11 
Lemma 2.2 Let f e R[x]. Then V{{x + 1)/) ^ V{f). 
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Proof : We proceed by induction in the number t of non-zero terms of /. The theorem is trivial 
for t = and t = 1. Now let us suppose that it holds for all t ^ n. Let / € M.[x] with 71 + 1 non-zero 
monomials. 



n+l 



/ — ttix"^ where 7^ for all i and ^ ai < 012 < • ■ • < 



deg(/) 



i=i 



Let g = Y^"=i CLix"' ■ By inductive hypothesis we have V{(x + l)g) ^ V{g). First, we consider 
the case a„ < d — 1, i.e., when the terms of {x + 1)5 do not overlap with those of a„+ix'*(x -I- 1). 
There are two posibilities: if a„a„+i > 0, then V{{x + 1)/) = V{{x + l)g) < V{g) = V{f), and 
if a„a„+i < 0, then V{{x + 1)/) = V{{x + l)g) + 1 ^ V{g) + 1 = V{f). In both cases we have 
V{{x + 1)/) ^ V{f)- Now it only remains the case a„ = d — 1. Here {x + 1)/ and (a; + l)g only 
differ in their terms of degree d and d -I- 1, as shown in the following table. 









(x + l)g 


an 





(^ + 1)/ 


O-n + 0,71+1 





If a„a„+i > 0, then V{f) ~ V{g), and according to the table, we have V{{x + 1)/) = V{[x + 1)5). 
Therefore V{{x + 1)/) < V{f). On the other hand, if a„a„+i < 0, then V{f ) = V{g) + 1, but we 
have three different posibilities for the table, depending whether |a„| is greater, equal or less than 
a„+i|. Set s = sgn(a„). 





x'^ 




{x + l)g 


s 







s 


—s 





x"^ 




(x + l)g 


s 










—s 





x"^ 




{x + l)g 


s 





{x + l)f 


— s 


—s 



Case |a„| > |a„+i 



Case |a„| = |a„+i| 



Case |a„| < |a„+i| 



The tables above show that V{{x + 1)/) < V{{x -I- \)g) + 1 for each of the three cases. Using the 
inductive hypothesis and V{f) — V{g) -f 1, we conclude that V{{x + 1)/) ^ V{f). □ 

Remark 2.3 Let f,gE M.[x] and suppose that g has t terms. Then V{f + g) ^ ^(/) + 2i. 

Note that the value of V^(0) is not relevant for thcorcm l2.2l The only reason for setting V{0) = —2 
is the previous remark. 

Proposition 2.4 Let f G R[x, y] be a polynomial with t non-zero terms. Then 

V{fix,x + 1)) sC 2t-2. 

Proof : We write / = J2^=i 0'i{x)y°"' , where < ai < • • • < a„ and ai{x) S R[x], and we set > 
the number of non-zero terms of a^. It is clear that t = ii + . . . + t„. 

We define fk = J27=k cii{x)y°''~°"' for fc = 1, . . . , n and /„+i — 0. Lemma [2?2] and Remark [2?3l imply 
that the polynomials fk satisfy: 







Vifkix,x + l)) ^ Vifk+i{x,x + l)) + 2tk 
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•.f = y"'fi =^ f{x,x+l) = {x+ir^h{x,x+l) =^ V{f{x,x+l))^V{h{x,x+l)). 
Thus, we conclude that V{f{x, x + 1)) ^ -2 + 2(ti + . . . + t„) = 2t - 2. □ 

Before finishing the proof of Theorcm ll.il we recaU Descartes' Rule of Signs: 

Theorem 2.5 (Descartes' rule of signs) Let f G M.[x] be a non-zero polynomial. Then f has at 
most V{f) positive roots counted with multiplicities. 

Proof of Theorem II .It If a = or 6 0, then g G M.[x] is a polynomial with at most t non-zero 
terms. Descartes' rule of signs implies that, either g = ot g has at most 2i — 1 ^ 6t — 4 real 
roots (counted with multiplicities except for the possible root 0). In the case a ^ and & ^ 0, 
the real roots of f{x, ax + h) correspond one to one to the roots of f{bx/a, b{x + 1)) = f{x, x + 1), 
where / = OiO""' a;"' j/'^'. Since this bijection preserves the multiplicity of the roots 

and maps the possible roots and ~b/a of g to the roots and —1 of /(a;, a; + 1), it sufhcies to 
consider the case a = b = 1, i.e. g = f{x,x + 1). Suppose that g ^ 0. Descartes' rule of signs 
and proposition 12.41 imply that the number of positive roots of g counting with multiplicities is at 
most 2t — 2. On the other hand, the roots of g in (— oo, —1) correspond to the positive roots of 
^ gi-l-x) = f{-l-x,-x) = fi{x,x + l), where /i = J2Liati-l)"'+^'x'^'y"''. Therefore the 
number of roots (with multiplicities) of g in {—oo, —1) is also bounded by 2t — 2. Finally, the roots 
of g in (—1,0) correspond to the positives roots of 

^ (. + 1)-^(^), (-^) . (. + 1)-^(^)/ {^,,^)- H^, ^ + 1) 

where /2 ~ J2i=i ai(— l)"'a;"'y'^°^'^^-'^"*^'^* . Therefore there are at most 2t — 2 of such roots. Taking 
into account the possible roots y — 1, counted each one at most once, we conclude that g has at 
most 6t — 4 real roots. □ 



3 Checking linear factors of a bivariate polynomial 

Proposition 3.1 Let f ~ J2l=i'^i^°'^y^^ ^ Let a, 6 G M such that \b\ / |1 — a|. Then 

y ~ ax — b \ f <^ — ax — 6 I f{x, a:") for at least 6i — 3 odd integers n ^ 3. 

Proof : (<;=) : Let 3 ^ ni < ?i2 < • • • < ngt-a the 6t — 3 odd numbers for which a;" — ax — 6 | f{x, x"). 
Let Wi GRhe a, root of x"* — aa; — 6 for each 1 ^ i ^ 6t — 3. Then f{wi, awi + b) = f{wi, w"*) — for 
all 1 ^ i < 6t — 3. This means that f{x, ax + b) has at least 6t — 3 real roots. Applying theorem [TTT] 
we conclude that f{x, ax + b) = 0, or simply y — ax — b \ f. It only remains to proof that Wi ^ wj 
for all i 7^ j. Actually, if a;"' — ax — b and x"^ ~ ax — b had a common root w = Wi = Wj G M, then 
^Ui-rij _ 2 g^j^^i therefore w = ±1. This would imply that — w"' — aw — b ^ — & ± (1 — a), in 
contradiction with the hypothesis |6| 7^ |1 — a|. □ 

Corollary 3.2 // / G M[a;, y] has t non-zero terms, then there is and odd integer 3^7i^ 12t — 5 
such that f{x, a:") ^ 0. 
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Proof : Otherwise every polynomial y — ax — b E M.[x,y] with |6| 7^ |1 — a| would divide /. □ 
Note that if (a, 6) ^ (0,±1), then either |6| 7^ |1 - a| or \b\ ^ \1 + a\. 

Algorithm TEST 

Input : A sparse polynomial / ~ J2l=i ciiX^'y^' G K[x,y] with t monomials, encoded as a list of 
vectors (a^, ai, Pi) £ A' x No x No representing the monomials of /, and two numbers a,b £ K. 

Output : True or False depending whether y — ax — b \ f{x, y) or not. 

Step 1 : If (a, 6) = (0, ±1), compute f{x,b). If this polynomial is zero, return True. Otherwise 
return False. 

Step 2 : If |fe| = |1 — a| then replace / by f{—x, y) and a by —a. 
Step 3 : For n = 3, 5, 7, . . . , 12t - 5 do 
Step 3.1 : If f{x, a;") ^ then 

Step 3.1.1 : Compute all the irreducible factors (with multiplicities) of x" — ax — in K[x] using 
a univariate dense factorization algorithm. 

Step 3.1.2 : Compute all the irreducible factors (with multiplicities) of f{x,x") in K[x] with 
degree n using a univariate sparse factorization algorithm. 

Step 3.1.3 : If there is an irreducible factor in the first list that either does not belong to the 
second list or belongs but with less multiplicity, then return False. 

Step 4 : Return True. 

The correctness of the algorithm is a consequence of the previous results. In order to estimate its 
complexity, we first state the following two famous results on the factorization of polynomials of 
univariate polynomials. 

DenseFactor Given / G K[x] of degree d and absolute height H, it is posible to compute all its 
irreducible factors in K[x] with multiplicities in [d[K : Q] logH]'^'-^^ bit operations (see 
for the rational case and [2] for the general case). 

SparseFactor Given f £ K[x] a. sparse polynomial of degree d, with at most t monomials and 
absolute height H, it is posible to find all its irreducible factors (with multiplicities) in K[x] 
of degree bounded by s in [ts[K : Q] \ogd logTJ]'^^^-' bit operations (see [Ij). 

The complexity of algorithm TEST is clearly dominated by its main loop (step 3), where it performs 
6i— 3 calls to DenseFactor and SparseFactor to factorize x^—ax—b completely and find all the factors 
of degree bounded by n of f{x, x"). We have that deg{x"—ax—b) = n ^ 12i— 5 and H{x^^—ax—b) ^ 
H{a)H{b), therefore step 3.1.1 requieres at most [{&t - 3) {Ut - 5) [K : Q] \og{H {a) H {b))f"^^^ bit 
operations. On the other hand, we have that /(x, x") is a sparse polynomial with at most t non-zero 
terms, of degree bounded by nd ^ {12t ~ 5)d and absolute height bounded by (2iJ)* because the 
coefficients of /(x,x") are sums of at most t coefficients of /. Thus, step 3.1.2 requieres no more 
than [{6t - 3)t{m - 5)[A : Q] log(d(12t - 5)) log(2i?)*]0(i) bit operations. This proves that the 
total number of bit operations performed by the algorithm is polynomial in t, \og{d), [K : Q] and 
\ogiHia)Hib)H). 
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